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Abstract 



We present an exact derivation of a process in which a microscopic 
measured system interacts with heat-bath and pointer modes of a measur- 
ing device, via a coupling involving a general Hermitian operator A of the 
system. In the limit of strong interaction with these modes, over a small 
time interval, we derive the exact effective many-body density matrix of the 
measured system plus pointer. We then discuss the interpretation of the 
dynamics considered as the first stage in the process of quantum measure- 
ment, eventually involving the wave-function collapse due to interactions 
with "the rest of the universe." We establish that the effective density ma- 
trix represents the required framework for the measured system and the 
pointer part of the measuring device to evolve into a statistical mixture de- 
scribed by direct-product states such that the system is in each eigenstate 
of A with the correct quantum-mechanical probability, whereas the expec- 
tation values of pointer-space operators retain amplified information of the 
system's eigenstate. 
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The problem of quantum measurement has fascinated scientists for a 
long time [1,2]. It has been argued that a large "bath" is an essential 
ingredient of the measurement process. Interaction with the bath, which 
might be a heat-bath in thermal equilibrium, causes decoherence which is 
needed to form a statistical mixture of eigenstates out of the initially fully 
or partially coherent quantum state of the measured system. An "exter- 
nal" bath ("the rest of the universe") may also play a role in selection of 
those quantum states of the pointer that manifest themselves in classical 
observations [2-7]. In this work we propose a model in which the pointer 
retains information on the measurement result because of its coupling to 
the measured system, without the need to couple it also to the internal 
bath. The measured system is still coupled to the internal bath. 

In an exactly solvable model of a quantum oscillator coupled to a heat 
bath of oscillators, it has been shown [4] that the reduced density matrix 
of the system, with the bath traced over, decoheres, i.e., it looses its off- 
diagonal elements in the eigenbasis of the interaction Hamiltonian. Recent 
work on decoherence [8-11] has explored the latter effect for rather general 
cases, for bosonic (oscillator) and spin baths. Applications for various phys- 
ical systems have been reported [12-18]. Fermionic heat bath has also been 
used in the literature [19]. 

It is clear, however, that the full function of a large, multimode mea- 
suring device, interacting with a small (microscopic) quantum system, must 
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be different from thermal equilibration or similar averaging effect. The de- 
vice must store and amplify the measurement outcome information. In this 
work we propose a solvable model that shows how this is accomplished. 

It must be stressed that for a complete description of the measurement 
process one needs to interpret the transfer of the information stored after 
the system-pointer and system-internal bath interaction to the macroscopic 
level [2]. Our attention here is no the process which corresponds to the 
first stage of the measurement, in which the pointer acquires amplified 
information by entanglement with the state of the system. Thus we do not 
claim to resolve the foundation-of-quantum mechanics issue of how that 
information is passed on to the classical world, involving the collapse of the 
wave functions of the system and each pointer mode. Indeed, it is unlikely 
that the wave function collapse can be fully described within the quantum- 
mechanical description of the three subsystems involved. Presumably, it 
would require consideration of an external bath with which the pointer and 
the internal bath interact. This problem is not presently solved [1-3], and we 
first sidestep it by assuming separation of time scales (see below). However, 
we later argue that our results provide useful hints on how to view the larger 
problem of quantum measurement. 

We now identify the three quantum systems involved. First, the mea- 
sured system, S, is a microscopic system with the Hamiltonian which will 
be also denoted by S. Second, the measuring device must have the "bath" 

-4- 



or "body" part, B, containing many individual modes. The kth mode will 
have the Hamiltonian B^. The bath part of the device is not observed, i.e., 
it can be traced over. Finally, the device must also have modes that are 
not traced over. These modes constitute the pointer, P, that amplifies the 
information obtained in the measurement process and can later pass it on 
for further amplification or directly to macroscopic (classical) systems. The 
rath pointer mode has the Hamiltonian P m . It is expected that expectation 
values of some quantities in the pointer undergo a large change during the 
measurement process. 

It turns out, a posteriori, that the device modes involved in the mea- 
surement process can be quite simple, and they need not interact with each 
other. This assumption allows us to focus on the evolution of the system 
S and its effect on the pointer P. However, it is the pointer's interaction 
with the external bath (some external modes, "the rest of the universe") 
that is presumed to select those quantum states of P that manifest them- 
selves classically. For now, let us avoid the discussion of this matter, see 
[2-6], by assuming that the added evolution of the pointer due to such ex- 
ternal interactions occurs on time scales larger than the measurement time, 
t. Similarly, when we state that the internal bath modes can be "traced 
over," we really mean that their interactions with the rest of the universe 
are such that these modes play no role in the wave-function-collapse stage 
of the measurement process. 
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Furthermore, the measurement process probes the wavefunction of the 
measured system at the initial time, t = 0, rather than its time evolution 
under S alone. It is ideally instantaneous. In practice, it is faster than the 
time scales associated with the dynamics under S. Such a process can be 
obtained as the limit of a system in which very strong interactions between S 
and B, and also between S and P, are switched on at t = and switched off 
at t > 0, with small time interval t. At later times, the pointer can interact 
with other, external systems to pass on the result of the measurement. 

Thus, we assume that the Hamiltonian of the system itself, S, can be 
ignored in the process. The total Hamiltonian of the system plus device 
will be taken as 

H = Y,Bk + Y, P ™ + bA J2 Lk+ P A 12 C ™ (!) 

km k m 

Here A is some Hermitian operator of the system that couples to certain 
operators of the modes, and C m . The parameters b and p are intro- 
duced to measure the coupling strength for the bath and pointer modes, 
respectively. They are assumed very large; the ideal measurement process 
corresponds to b,p — > oo. 

We note that the modes of P and B can be similar. The only difference 
between the bath and pointer modes is in how they interact with the "rest 
of the universe": the bath is traced over (unobserved), whereas the pointer 
modes have their wave functions collapsed in a later step of the measurement 
process. Thus, we actually took the same coupling operator A for the bath 
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and pointer. In fact, all the exact calculations reported in this work can 
be also carried out for different coupling operators A& and A p , for the bath 
and pointer modes, provided they commute, [A&, A p ] = 0, so that they share 
a common set of eigenfunctions. The final wavefunction of the measured 
system, after the measurement, is in this set. Analytical calculation can be 
even extended to the case when the system's Hamiltonian S is retained in 
(1), provided all three operators, S, A&, A p , commute pairwise. The essential 
physical ingredients of the model are captured by the simpler choice (1). 

We will later specify all the operators in (1) as the modes of the bosonic 
heat bath of Caldeira-Leggett type [17,19-26]. For now, however, let us 
keep our discussion general. We will assume that the system operator A 
has nondegenerate, discrete spectrum of eigenstates: 

A|A) = A| A) (2) 

Some additional assumptions on the spectrum of A and S will be encoun- 
tered later. We also note that the requirement that the coupling parameters 
b and p are large may in practice be satisfied because, at the time of the 
measurement, the system's Hamiltonian S corresponds to slow or trivial 
dynamics. 

Initially, at t — 0, the quantum systems (S, B, P) and their modes are 
not correlated with each other. We assume that p is the initial density ma- 
trix of the measured system. The initial state of each bath and pointer mode 
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will be assumed thermalized, with (3=1/ (kT) and the density matrices 



,-0 B k 



Ob = 



Tr fc (e-^*) 



Tr m (e-Pr™) 



(3) 



We cannot offer any fundamental physical reason for having the initial bath 
and pointer mode states thermalized, especially for the pointer; this choice 
is made to allow exact solvability. 

The density matrix of the system at the time t is 



R = e- im ' n p(l[e h ) (j[a, 

V k ) \ m 



iHt/h 



(4) 



The bath is not probed and it can be traced over. The resulting reduced 
density matrix r of the combined system S + P will be represented by 
its matrix elements in the eigenbasis of A. These quantities are each an 
operator in the space of P: 



r X y = (X\Tt b (R)\X') 



(5) 



We now assume that operators in different spaces and of different modes 
commute. Then one can show that 



r\\> = p\\> 



~Y^ e - %t ^ p ^^c m )lh G ^ e it{p m ^ v \'c m )lh 



m 



X 



Y[ Tr fc j^-^Bk+bXLk) /hQ ke it(B k +b\' L k ) /h 

. k 



(6) 



where p\\> = (\\p\X'). This result involves products of P-space operators 
and traces over B-space operators which are all single-mode. Therefore, 
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analytical calculations are possible for some choices of the Hamiltonian (1). 
The observable A can be kept general. 

The role of the product of traces over the modes of the bath in (6) is 
to induce decoherence which is recognized as essential for the measurement 
process, e.g., [1,2]. At the time t, the absolute value of this product should 
approach 5\\> in the limit of large b. Let us now assume that the bath 
is bosonic. The Hamiltonian of each mode is then TiUka\ak, where for 
simplicity we shifted the zero of the oscillator energy to the ground state. 
The coupling operator is usually selected as Lk = gl a k + 9k%- For 
simplicity, though, we will assume that the coefficients g\~ are real: 



For example, for radiation field in a unit volume, coupled to an atom [27], 
the coupling is via a linear combination of the operators {a^ + cl\)/ y/uJk and 



and momentum, respectively. Our calculations can be extended to have an 
imaginary part of which adds interaction with momentum. 

The product of traces in (6) can be calculated by coherent-state or 
operator- identity techniques [8-10]. Here and below we only list the results 
of such calculations which are usually quite cumbersome: 



J] Tr fc {...} = exp {-2b 2 (A - A') 2 T(t) + ib 2 [A 2 - (A') 2 ] 7 (*)} (8) 






i{ak — a\)/ y/uJk- For a spatial oscillator, these are proportional to position 



k 




(9) 



k 
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Explicit form of 7(f) is also known [8]. 

In the continuum limit of many modes, the density of the bosonic bath 
states in unit volume, D(u>), and the Debye cutoff with frequency, ujd, are 
introduced [22] to get 








Let us consider the popular choice termed Ohmic dissipation [22] , motivated 
by atomic-physics [27] and solid-state applications [22], corresponding to 

D(co)g 2 (u) = ttu (11) 

where Q is a constant. Other powers of u have also been considered, e.g., 
[11]. In studies of decoherence [8-11] for large times t, for models without 
strong coupling, not all the choices of D(u)g 2 (u) lead to complete decoher- 
ence [11] because T(t) must actually diverge to +00 for t ^> h0, as happens 
for the choice (11). 

Let us assume that the energy gaps of S are bounded so that there 
exists a well defined time scale h/ AS of the evolution of the system under S. 
There is also the time scale 1/ljd set by the frequency cutoff assumed for the 
interactions. The thermal time scale is Ti/3. The only real limitation on the 
duration of measurement is that t must be less then h/ AS. In applications, 
typically [22] one can assume that 1/cjd <C h/AS. Furthermore, it is 
customary to assume that the temperature is low [22], 

t and l/u D <C Ti/AS <C Tif3 (12) 



In the limit of large h(3, the absolute value of (8) reduces to 



Abs|nTr fc {...}j ^expj-A^^-AO'Ml + ^t) 2 ]! (13) 

In order to achieve effective decoherence, the product (A A) 2 b 2 lnfl + ^^t) 2 ] 
must be large. The present approach only applies to operators A with 
nonzero scale of the smallest spectral gaps, AA. 

We note that the decoherence property needed for the measurement 
process will be obtained for nearly any well-behaved choice of Dfa)g 2 fa) 
because we can rely on the value of b being large rather than on the proper- 
ties of the function T(t). If b can be large enough, very short measurement 
times are possible. However, it may be advisable to use measurement times 
1/Vd <t< h/ AS to get the extra amplification factor ~ Infant) and al- 
low for fuller decoherence and less sensitivity to the value of t in the pointer 
part of the dynamics, to be addressed shortly. We notice, furthermore, 
that the assumption of a large number of modes is important for monotonic 
decay of the absolute value of (8) in decoherence studies [8-11], where ir- 
reversibility is obtained only in the limit of infinite number of modes. In 
our case, it can be shown that such a continuum limit allows to extend the 
possible measurement times from £<C 1/Vd to 1/Vd < h/ AS. 

Consider the reduced density matrix r of S + P, see (6). It becomes 
diagonal in |A), at the time t, because all the nondiagonal elements are 



small, 

r = ^ |A)(A| pxx Y[e- u ^ +pXC ^ /h a m e u ^ +pXC ^ /h (14) 

A m 

Thus, the described stage of the measurement process yields the density ma- 
trix that can be interpreted as describing a statistically distributed system, 
without quantum correlations. This, however, is only meaningful within 
the ensemble interpretation of quantum mechanics. 

For a single system plus device, coupling to the rest of the universe is 
presumably needed (this problem is not fully understood in our opinion, see 
[2]) for that system to be left in one of the eigenstates |A), with probability 
px\. After the measurement interaction is switched off at t, the pointer 
coupled to that system will carry information on the value of A. This 
information is "amplified," owing to the large parameter p in the interaction. 

We note that one of the roles of the pointer having many modes, many 
of which can be identical and noninteracting, is to allow it (the pointer 
only) to still be treated in the ensemble, density matrix description, even if 
we focus on the later stages of the measurement when the wave functions of 
a single measured system and of each pointer mode are already collapsed. 
This pointer density matrix can be read off (14). This aspect is new and 
it may provide a useful hint on how to set up the treatment of the full 
quantum-measurement process description. 

Another such hint is provided by the fact that, as will be shown shortly, 
the changes in the expectation values of some observables of the pointer re- 



tain amplified information on the system's eigenstate. So, coupling to the 
rest of the universe that leads to the completion of the measurement process, 
should involve such an observable of the pointer. Eventually, the informa- 
tion in the pointer, perhaps after several steps of amplification, should be 
available for probe by interactions with classical devices. 

At time t = 0, expectation values of various operators of the pointer 
will have their initial values. These values will be different at the time t 
of the measurement owing to the interaction with the measured system. It 
is expected that the large coupling parameter p will yield large changes in 
expectation values of the pointer quantities. This does not apply equally 
to all operators in the P-space. Let us begin with the simplest choice: the 
Hamiltonian ^ P m of the pointer. 

m 

We will assume that the pointer is described by the bosonic heat bath 
and, for simplicity, use the same notation for the pointer modes as that used 
for the bath modes. The assumption that the pointer modes are initially 
thermalized, see (3), was not used thus far. While it allows exact analytical 
calculations, it is not essential: the effective density matrix describing the 
pointer modes at the time t, for the system's state A, will retain amplified 
information on the value of A for general initial states of the pointer. 

This effective density matrix is the product over the P-modes in (14). 
For the "thermal" <j m from (3), the expectation value of the pointer energy 



Ep can be calculated from 



(Ep)\ Tip ( e -** = Tr P | ^^ m a^a m j x (15) 



n 



The right-hand side can be reduced to calculations for individual modes. 
Operator identities can be then utilized to obtain the results 

(E P ) x (t) = (E P )(0) + (AEp)x(t) (16) 

(E P )(0) = fiJ2^me- hpuJm (1 - e - hf3 ^y 2 (17) 



m 



For a model with Ohmic dissipation, the resulting integral, in the continuum 
limit, can be calculated to yield 

<«*«> - m r v . ;?:;■ 

which should be compared to the exponent in (13). The energy will be 
an indicator of the amplified value of the square of A, provided p is large. 
Furthermore, we see here the advantage of larger measurement times, t ^> 
1/lod- The change in the energy then reaches saturation. After the time 
t, when the interaction is switched off, the energy of the pointer will be 
conserved. 
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Let us consider the expectation value of the following Hermitian oper- 
ator of the pointer: 

m m 

For an atom in a field, X is related to the electromagnetic field operators 
[24]. One can show that (X P )(0) = and 



(AX P ) x (t) = (X P )x(t) = 




~ h 1 + (to D t) 2 

The change in the expectation value of X is linear in A. However, this 
operator is not conserved. One can show that after the time t its expectation 
value decays to zero for times t + 0{\/ujd)- 

We note that by referring to "unit volume" we have avoided the dis- 
cussion of the "extensivity" of various quantities. For example, the initial 
energy (Ep)(0) is obviously proportional to the system volume, V. How- 
ever, the change (AEp)\(i) will not be extensive; typically, g 2 {oj) oc 1/V, 
D(uj) oc V. Thus, while the amplification in our measurement process 
can involve a numerically large factor, the changes in the quantities of the 
pointer will be multiples of microscopic values. Multi-stage amplification, 
or huge coupling parameter p, would be needed for the information in the 
pointer to become truly "extensive" macroscopically. 



In practice, there will be probably two types of pointer involved in a 
multistage measurement process. Some pointers will consist of many nonin- 
teracting modes. These pointers carry the information, stored in a density 
matrix rather than a wave function of a single system. The latter transfer- 
ence hopefully makes the wavefunction collapse and transfer of the stored 
information to the macroscopic level less "mysterious and traumatic." The 
second type of pointer will involve strongly interacting modes and play the 
role of an amplifier by utilizing the many-body collective behavior of the 
coupled modes (phase-transition style). Its role will be to alleviate the artifi- 
cial requirement for large mode-to-system coupling parameters encountered 
in our model. 

In summary, we described the first stage of a measurement process. It 
involves decoherence due to a bath and transfer of information to a large 
system (pointer) via strong interaction over a short period of time. The 
pointer itself need not be coupled to the internal bath. While we do not offer 
a solution to the foundation-of-quantum-mechanics wave-function collapse 
problem [2], our results do provide two interesting observations. 

Firstly, the pointer operator "probed" by the rest of the universe during 
the wave-function collapse stage, may be in part determined not only by 
how the pointer modes are coupled to the external bath [3-7], but also by 
the amplification capacity of that operator in the first stage of the process, 
as illustrated by our calculations. 

- 16 - 



Secondly, for a single system (rather then an ensemble), the multiplicity 
of the (noninteracting) pointer modes might allow the pointer to be treated 
within the density matrix formalism even after the system and each pointer- 
mode wave functions were collapsed. Since it is the information in the 
pointer that is passed on, this observation might seem to resolve part of the 
measurement puzzle. Specifically, it might suggest why only those density 
matrices entering (14) are selected for the pointer: they carry classical 
(large, different from other values) information in expectation values, rather 
than quantum-mechanical superposition. However, presumably [2] only a 
full description of the interaction of the external world with the system 
S + P can explain the wavefunction collapse of S. 

We acknowledge helpful discussions with Professor L. S. Schulman. 
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grant DAAD 19-99-1-0342. 
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